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Abstract 

A study of W + W~ events accompanied by hard photon radiation, E 1 > 2.5 GeV, produced in 
e + e~ collisions at LEP is presented. Events consistent with being two on-shell W-bosons and an 
isolated photon are selected from 681 pb _1 of data recorded at 180 GeV < <Js < 209 GeV. From 
the sample of 187 selected W + W~7 candidates with photon energies greater than 2.5 GeV, the 
W + W~7 cross-section is determined at five values of ^fs. The results are consistent with the 
Standard Model expectation. Averaging over all energies, the ratio of the observed cross-section 
to the Standard Model expectation is 

i?(data/SM) = 0.99 ± 0.09 ± 0.04, 

where the errors represent the statistical and systematic uncertainties respectively. These data 
provide constraints on the related 0(a) systematic uncertainties on the measurement of the 
W-boson mass at LEP. Finally, the data are used to derive 95 % confidence level upper limits 
on possible anomalous contributions to the W + W~77 and W + W~Z°7 vertices: 

-0.020 GeV" 2 < a /A 2 < 0.020 GeV" 2 , 
-0.053 GeV" 2 < a c /A 2 < 0.037 GeV" 2 , 
-0.16 GeV" 2 < a n /A 2 < 0.15 GeV" 2 , 

where A represents the energy scale for new physics and do, a c and a n are dimensionless coupling 
constants. 
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1 Introduction 

The W + W~ pair production cross-section has been precisely measured at LEP over a range 
of centre-of-mass energies [1-3]. The data are well described by the Standard Model (SM) 
expectation [4, 5] . The good agreement between experiment and theory is only obtained once 
factorizable and non-factorizable 0(a) photonic corrections are included in the theoretical cal- 
culations (see for example [6] and references therein). The inclusion of real and virtual photonic 
corrections in the YFSWW [4] and RacoonWW [5] programs has reduced the theoretical un- 
certainty on the CC03 e + e~ — > W + W~ cross-section 1 to below 0.5 % [6]. Uncertainties in these 
0(a) corrections may lead to small, but non-negligible, systematic uncertainties in the deter- 
mination of the W-boson mass, M w , at LEP [7]. This paper presents a study of the process 
e + e~ — > W + W~7 and thus probes the modelling of real photonic corrections to the W + W~ 
pair creation process. The data are used to obtain measurements of the e + e~ — > W + W~7 
cross-section within a restricted phase-space region, <xww7> f° r 180 GeV < yfs < 209 GeV. 

In the SM, photon radiation in the W + W~ production process at LEP can be categorized 
into four main classes of diagrams: initial state radiation (ISR); final state radiation (FSR) 
from a lepton; FSR from the quark or from the associated parton shower; and bremsstrahlung 
from one of the intermediate W-bosons, referred to as WSR. At LEP energies WSR has a sig- 
nificant effect only through interference with ISR. Experimentally photons arising from decays 
of hadrons in a jet are indistinguishable from FSR photons from a quark or parton shower. 
For this reason, and due to the relatively large uncertainties in the Monte Carlo modelling of 
photon production in the parton shower, all photons associated with hadronic jets (from hadron 
decay and FSR) are considered background for the measurements of <Jww7- 

The measurements of the W + W~7 cross-section are compared with the predictions of the 
KORALW [8], KandY [9] (the concurrent Monte Carlo KORALW1.51 and YFSWW3) and 
RacoonWW [5] programs. These comparisons are used to obtain the first data-driven estimate 
of the systematic uncertainty on M w due to the Monte Carlo description of real photon radiation 
in W + W~ events. 

In addition, the W + W~7 final state is sensitive to possible anomalous W + W~77 and 
W + W~Z°7 quartic gauge boson couplings (QGCs). At LEP energies the contribution of the 

1 CC03 refers to the three doubly resonant diagrams for e + e~ — > W + W~. 
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SM QGC diagram is negligible. The data presented in this paper are used to place upper 
limits on the size of possible anomalous QGCs. These limits are more than a factor three 
tighter than previous OPAL results from e + e~ — > W + W~7 [10] and are consistent with other 
measurements [11]. 

2 The OPAL Detector, Data Samples and Monte Carlo 

2.1 The OPAL Detector 

The OPAL detector includes a 3.7 m diameter tracking volume within a 0.435 T axial mag- 
netic field. The tracking detectors include a silicon micro-vertex detector, a high precision gas 
vertex detector and a large volume gas jet chamber. The tracking acceptance corresponds to 
approximately | cos#| < 0.95 (for the track quality cuts used in this study) 2 . Lying outside the 
solenoid, the electromagnetic calorimeter (ECAL) consisting of 11 704 lead glass blocks has full 
acceptance in the range | cos^l < 0.98 and a relative energy resolution of approximately 6 % 
for 10 GeV photons. The magnet return yoke is instrumented with streamer tubes which serve 
as the hadronic calorimeter. Muon chambers outside the hadronic calorimeter provide muon 
identification in the range | cos#| < 0.98. A detailed description of the OPAL detector can be 
found in [12]. 

2.2 Data Sample 

During LEP2 operation the centre-of-mass energy was increased from 161 GeV to 209 GeV 
in several steps. The total integrated luminosity of the data sample considered in this paper, 
evaluated using small angle Bhabha scattering events observed in the silicon tungsten forward 
calorimeter [13], is (681 ± 2) pb _1 . For the purpose of measuring the W + W~7 cross-section 
these data are divided into the five yfs ranges listed in Table 1. These ranges reflect the main 
energy steps as the centre-of-mass energy was increased during LEP2 operation. The data 
recorded at 161 GeV and 172 GeV, corresponding to a total integrated luminosity of 20 pb~ x , 
are not used here. 



Range/GeV 


< y/s >/GeV 


r/pb- 1 


180.0-185.0 


182.68 


57.2 


188.0-189.0 


188.63 


183.1 


191.0-196.0 


194.44 


105.7 


199.0-204.0 


200.21 


114.1 


204.0-209.0 


205.92 


220.6 



Table 1: The energy binning used for the W + W~7 cross-section measurements. The y^s range 
covered by each bin, the mean luminosity weighted value of y/s and the corresponding integrated 
luminosity, C, are listed. 



2 The OPAL right-handed coordinate system is defined such that the origin is at the centre of the detector 
and the z axis points along the direction of the e~ beam; 8 is the polar angle with respect to the z axis. 
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2.3 Monte Carlo 



A number of Monte Carlo (MC) samples, all including a full simulation [14] of the OPAL 
detector, are used to simulate the SM signal and background processes. For this paper the main 
MC samples for the process e + e~ — > W + W~7 were generated using the KandY [9] program and, 
unless otherwise specified, the SM expectations for the e + e~ — > W + W~7 cross-section refer to 
the KandY prediction. KandY includes exact O(a) YFS exponentiation [15] for the W + W~ 
production process, with 0(a) electroweak non-leading (NL) corrections combined with YFS 
exponentiated 0(a 3 ) leading logarithm (LL) initial state radiation. Final state radiation from 
leptons is implemented in PHOTOS [16] and radiation from the quark induced parton-shower 
is performed by JETSET [17]. The most notable improvements over the KORALW program 
are the leading non-factorizable corrections in the Screened Coulomb ansatz [18], the inclusion 
of bremsstrahlung from the W-pairs (WSR), and the implementation of O(a) electroweak NL 
corrections. 

The KORALW program [8] is used to simulate the background from four-fermion final states 
which are incompatible with coming from the decays of two W-bosons (e.g. e + e~ — > qqu. + /U~). 
The two-fermion background processes, e + e~ — > Z°/7 — > qq and e + e~ — > Z°/7 — > t + t~, 
are simulated using KK2F [19]. The background in the W + W~7 event selection from multi- 
peripheral two-photon diagrams was found to be negligible. 

In addition, the RacoonWW program [5] is used in the Improved Born Approximation 
(IBA) mode to obtain independent predictions of the cross-sections for e + e~ — > W + W~7 and 
e + e~ — > 4/7. In this mode all lowest order diagrams contributing to e + e~ — > W + W~7 are 
included. The EEWWG program [20] is used to obtain predicted cross-sections in the presence 
of anomalous QGCs which are then used to extract experimental limits on the anomalous 
contributions to the W + W~77 and W + W~Z°7 vertices. 

3 W+W~7 Signal Definition 

The process e + e~ — > W + W~7 results in a four-fermion plus photon final state, /1/2/3/47) 
where the fermion flavours are appropriate for W-decay. In the SM, photons are radiated in 
several classes of diagrams corresponding to ISR, FSR from both charged leptons and quarks, 
radiation from the W-boson (WSR) and the Standard Model QGC diagram. The invariant 
mass distributions of the fermions are different for the different radiation processes. In the case 
of ISR, the /i/ 2 and fsf^ systems are produced with invariant masses close to M\y. In the case 
of FSR, the /i/ 2 7 and / 3 / 4 combinations or the f 1 f 2 and fsf^y combinations give invariant 
masses close to the W-boson mass. For photon energies E 1 > T w , where T w is the W-boson 
width, events from FSR tend to occupy a different kinematic region from those arising from 
the ISR or QGC diagrams. Consequently, interference between FSR and ISR/QGC diagrams 
is suppressed. At LEP energies the effect of WSR is only significant through interference with 
ISR; the WSR diagrams are only of relevance to the region of phase-space populated by ISR 
diagrams. 

Only part of the W + W~7 phase-space is accessible experimentally and, therefore, it is 
necessary to define a specific region of phase-space in which the cross-section will be measured. 
The definition of the signal region is chosen to be well matched to the experimental sensitivity. 
In addition, by defining the cross-section to correspond to a region of four-fermion phase- 
space dominated by the doubly resonant W + W" production (CC03) diagrams, contributions 
from other interfering diagrams can be made small. In this way, the experimental results can be 
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compared with both the predictions of calculations implementing all four-fermion diagrams and 
with calculations implementing only CC03 diagrams. Finally, invariant mass cuts are imposed 
to reduce the contribution of FSR both from quarks and from leptons. This is desirable for 
two reasons. Firstly, any new physics is unlikely to manifest itself in a modification of FSR. 
Secondly, it reduces modelling uncertainties which are potentially large in the case of FSR from 
the quark-induced parton shower. 

In this paper, the W + W~7 — ► /1/2/3/47 cross-section, denoted by (JWW7, is measured for: 

• Ej > 2.5 GeV, where E 1 is the photon energy. 

• |cos# 7 | < 0.975, where cos6> 7 is the cosine of the polar angle of the photon. 

• cos < 0.90, where cos# 7 / is the cosine of the minimum angle between the photon and 
any of the charged fermions in the four-fermion final state. 

• |cos^| < 0.95, where |cos^| is the modulus of the cosine of the polar angle of the 
charged lepton in the W + W~ — > qqiu^ final state. In the W + W~ — > t^v^Vi final state 
this requirement applies to both of the charged leptons. 

• — Mw| and \M^ — Mw| < 3 IV, where M^j and A^/ 3 j 4 are the invariant masses 
of fermions consistent with being from the decays of the W~ or W + . 

The first three requirements are closely matched to the ability to reconstruct a pure sample 
of isolated photons in the OPAL detector. The requirement on the polar angle of the charged 
leptons from W-decay is imposed because the W + W~ event selection becomes significantly less 
efficient beyond the acceptance of the tracking chambers. It also reduces contributions from 
interfering four-fermion background diagrams such as the t-channel process e + e~ — > WeV e . The 
cut on the invariant masses of the fermion pairs further reduces the (interfering) four-fermion 
backgrounds and suppresses the contribution of FSR to the signal region. Due to the finite jet 
width, jets are detected over the full polar angle acceptance and therefore there is no explicit 
requirement on the polar angle of the quark. 

In the above definition of the signal, all requirements are made on generator level quantities. 
Generator level refers to the true four-momenta of particles in the /1/2/3/47 final state. The 
cross-section within the above kinematic cuts, <3"ww7, is dominated by doubly-resonant W + W~ 
production. For example, the difference between the cross-section for the full set of 4/7 dia- 
grams relative to cross-section for the CC03 diagrams alone is less than 0.5 % (calculated using 
the IB A implemented in RacoonWW [5]). 

4 W + W~7 Event Selection 

The selection of W + W~7 events proceeds in three stages: selection of W + W~ events, photon 
identification, and background rejection using kinematic information. All W + W~ final states 
are used in this study. 

4.1 W+W- Selection 

The W + W~— >£ + v e l~V£, W + W~ ^qqiV e and W + W _ ^qqqq selections of reference [1] are used 
as the basis of the W + W~7 selections 3 . For W + W~(7) -^i^v^'Vf^ and W + W~(7) — >qq^7 

3 Reference [1] refers to the event selection at y/s = 189 GeV. For data recorded at higher centre-of-mass 
energies the same likelihood selection is used but with reference distributions obtained from Monte Carlo events 
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the standard selections are applied. For W + W~(7) — >qqqq~7 events, a modified version of the 
W + W~— >qqqq selection of reference [1] is used. In the standard selection, events are forced 
into four jets using the Durham algorithm [21]. In approximately 10 % of Monte Carlo 
events with high energy photons {E 1 > 10 GeV), the photon alone forms one of the four jets. 
This introduces an additional inefficiency, due to the requirement in the preselection that there 
should be at least one charged particle track associated with each jet. For this reason, events 
failing the standard W + W~^qqqq selection are forced into four jets after excluding the highest 
energy isolated electromagnetic calorimeter cluster and the selection re-applied. The overall 
selection efficiency for W + W~7 events within the signal definition is 88 % and is approximately 
independent of centre-of-mass energies for 180 GeV < y/s < 209 GeV. 

4.2 Photon Identification 

Photon identification is similar to that described in [22], although for this study the minimum 
photon energy is reduced to 2.5 GeV. Photon candidates are identified as one of three types: 

• Unassociated ECAL clusters defined by the requirement that no charged particle track, 
when extrapolated to the front-face of the ECAL, lies within a distance defined by the 
typical angular resolution of the ECAL cluster. The lateral spread of the cluster was 
required to satisfy the criteria described in reference [22] . 

• Two-track photon conversions which are selected using an artificial neural network as 
described in [23]. 

• Conversions where only a single track is reconstructed, identified as an electromagnetic 
calorimeter cluster associated with a track which is consistent with originating from a 
photon conversion. The track is required to have no associated hits in either layer of the 
silicon micro-vertex detector or in the first six layers of the central vertex chamber. 

For both types of conversion, the photon energy is defined by the sum of cluster energies pointed 
to by the track(s). 

Photon candidates identified using the above criteria are required to satisfy isolation re- 
quirements. The summed energies of any additional tracks and clusters in a 20° half-angle cone 
defined by the photon direction have to be less than 2 GeV. In addition, the energy deposited 
in the hadron calorimeter in a 20° half-angle cone around the photon candidate is required to 
be less than 5 GeV. If the invariant mass formed from the photon candidate and the energy 
deposit in any ECAL cluster is less than 0.25 GeV/c 2 the candidate is rejected in order to 
suppress photons from n° decay. For photon candidates with 2.5 GeV < E y < 10.0 GeV a rel- 
ative likelihood selection is applied to reduce the background from photons from the decays of 
hadrons (dominated by 7r° and rj decays). The likelihood is based on five discriminant variables: 
Ej, | cos 9 1 1, the angle between the photon and the nearest jet, the angle between the photon 
and the nearest track, and the minimum invariant mass formed from the photon candidate and 
any other ECAL cluster in the event. For photons above 10 GeV the background is low and 
no photon identification likelihood is needed. 

4.3 Photon Acceptance 

The identified photon is required to lie within the polar acceptance, 
generated at higher centre-of-mass energies. 
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• |cos6» 7 | < 0.975. 

The photon is also required to be isolated from the charged fermions in the final state. Cuts are 
applied on the cosine of the angle between the photon and closest jet, cos # 7 _ jet, and on the 
cosine of the angle between the photon and a charged lepton from the W-boson decay, cos# 7 ^: 

• cos # 7 _ jet < 0.9 for W + W _ — >qq£z/£7 and W + W~^qqqqy events, 

• cosft^ < 0.9 for W + W~^£ + z/^~77£7 and W + W~^qq£z7^7 events. 

For selected events with photons within the generator level acceptance the photon identi- 
fication efficiency is 75 % for E 1 > 7.5 GeV, 69 % for 5.0 GeV < E~ f < 7.5 GeV and 45 % for 
2.5 GeV < E 1 < 5.0 GeV. The photon identification efficiency is almost independent of cos# 7 in 
the region |cos# 7 | < 0.975. The non-photonic backgrounds are less than 4 % for |cos6> 7 | < 0.95. 
For |cos# 7 | > 0.95 the background increases to 8 %. If more than one photon candidate passes 
the photon acceptance requirements only the highest energy photon is retained for the following 
analysis. 

4.4 Kinematic Requirements 

The photon in selected W + W~7 events is classified as ISR, FSR from the lepton, or as being 
associated with a jet (either FSR from the parton shower or coming from hadron decay). 
No special treatment is made for WSR because WSR diagrams are only observable through 
interference with ISR diagrams and, consequently, the effects of WSR diagrams will be apparent 
in the event sample classified as ISR. In £ + u e £~V e events, photons are classified as ISR if 
cos# 7 ^ < |cos0 7 |, otherwise the photons are classified as FSR from one of the charged leptons. 
For the qq^v and qqqq channels the classification is performed using a relative likelihood 
selection in which kinematic fitting plays a major role. Three kinematic fits are employed, 
corresponding to the following hypotheses: 

a) the photon originates from FSR from the quark; the fit assumes a two-body W + W~ final 
state, where the identified photon is included as part of the nearest jet. 

b) the photon originates from FSR from the lepton (only used for W + W~(7) -^qqiu^ 
events); the fit assumes a two-body W + W~ final state, where the photon is associated 
with the charged lepton. 

c) the photon originates from ISR; the fit assumes a three body final state consisting of the 
two W-bosons and the photon. 

In each case, the constraints of energy and momentum conservation are imposed and the two 
reconstructed masses of the W-boson candidates are required to be equal [24]. An event is 
considered consistent with one of the above hypotheses if the fit converges with a fit probability 
of greater than 0.1 % and if the reconstructed W-boson mass is greater than 74 GeV. In fully 
hadronic events there are three possible jet-pairing combinations. Here, for each fit hypothesis, 
the combination yielding the highest kinematic fit probability is used. 

The reconstructed W-boson mass from the three kinematic fit hypotheses along with the 
cosine of the angle between the photon and the nearest jet are used as the inputs to the relative 
likelihood. For qq£z7^ events the cosine of the angle between the photon and the charged 
lepton is also used. The distributions used in the relative likelihood classification are shown in 
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Figure 1. Good agreement between data and simulation is observed. Three relative likelihoods 
are constructed and events are classified as being either from ISR, FSR from the charged lepton 
or radiation associated with the jets. The resulting ISR relative likelihood distribution, Cisr 
is shown in Figure If. Events are classified as ISR if Cisr > Cfsr and Cisr > Cjet, where 
Cfsr and Cjet are the relative likelihoods for the respective hypotheses of FSR and radiation 
associated with the jets. Only those W + W~7 candidate events classified as ISR are retained 
for the analysis. These events are consistent with on-shell W-bosons (fit c), M^j ~ ^fzf A ~ 
Mw and an isolated photon. This procedure suppresses events with final state radiation and 
events where the photon is from hadron decay. It also significantly reduces background from 
e + e~ — > qqy. As a result the systematic uncertainties from photons associated with jets (FSR 
and tc°/t] decays) are greatly reduced. 

The application of the above kinematic requirements retains approximately 75 % of selected 
signal W + W~7 events with an identified photon (using the definition of Section 3) whilst 
rejecting 85 % — 98 % (increasing with the photon energy) of events with photons either from 
FSR or from the decays of mesons. 



5 Measurement of the W + W 7 cross-section 

Using the selection criteria defined in the previous section, 187 W + W~7 events with E 1 > 
2.5 GeV are selected compared to the KandY expectation of 188.4 ± 1.0 events (where the 
error on the expectation is the quadrature sum of the MC statistical error and luminosity error). 
Figure 2a shows the photon energy spectrum for the selected W + W~7 events. Figure 2b shows 
the distribution of |cos 7 | and Figure 2c shows the distribution of the cosine of the angle between 
the photon and the nearest charged fermion from the reconstructed W-decay (i.e. lepton or jet) 
in the event. Good agreement between data and Monte Carlo is observed for all distributions. 
The effect of an anomalous QGC on the photon energy and polar angle distributions is also 
shown. 

The W + W~7 cross-section is determined within the acceptance defined in Section 3 for the 
five mean centre-of-mass energies listed in Table 1. The W + W~7 cross-section is calculated 
from 

(Nobs - <t B gd£) 

<T WW7 — — , 

CWW7^WW7Z- 

where iV bs is the accepted number of events, <tbgd is the SM background cross-section and C is 
the integrated luminosity. The selection efficiency for events generated within the acceptance 
defined in Section 3, £\yw7, is evaluated using KandY MC W + W~7 events. Background from 
migration of W + W~7 events from just outside the signal region into the selected event sample 
due to finite detector resolution is accounted for by a factor cwwi- This allows the contribution 
from selected W + W~7 events outside the signal definition but within the acceptance E 1 > 
2.0 GeV and |cos # 7 | < 0.98 to scale with the measured cross-section (in contrast to treating this 
component as background which is fixed by Monte Carlo expectation). The selection efficiency, 
£ww7> varies from 41 % — 47 % increasing with centre-of-mass energy. The correction factor, 
C WW7> is 114 and is almost independent of centre-of-mass energy. The background cross- 
section, <Jbgd, is estimated using KandY and KK2F. The background from W + W~ events 
with photons associated with the jets, including photons from FSR from the parton-shower, is 
scaled by a factor of 1.30±0.15, as described in Section 5.1, to account for known discrepancies 
between data and the JETSET prediction. The Monte Carlo predicts that 24 % of the selected 
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event sample arises from background processes (including photons associated with jets). The 
results are listed in Table 2 where they are compared to the predictions from KandY, and 
are displayed in Figure 3. The systematic uncertainties are described in Section 5.1. For the 
purpose of combination with the other LEP experiments, the results for a more restrictive signal 
acceptance are given in the Appendix. 





0"WW7/fb> 


< v/i >/GeV 


Data 


KandY 


182.68 


277 ±117 ±13 


327 ±3 


188.63 


388 ± 74 ±17 


378 ±4 


194.44 


255 ± 84 ±15 


411 ±4 


200.21 


459 ±100 ±20 


427 ±4 


205.92 


489 ± 73 ±21 


443 ±4 



Table 2: W + W~7 cross-section measurements for the five centre-of-mass energies listed in Table 
1. The errors on the measurements are statistical and systematic respectively. The errors on 
the KandY expectations are due to limited Monte Carlo statistics. 

Table 3 shows the ratio of measured to predicted W + W~7 cross-sections averaged over 
the five values of y/s for the theoretical predictions from KandY, RacoonWW, EEWWG and 
KORALW. For RacoonWW, FSR from the parton shower is included as signal since, unlike 
for KandY and KORALW, there is no way of removing its contribution at the generator level. 
As a consequence of the uncertainties of the modelling of photons from the parton shower 
this results in an increased systematic uncertainty as discussed in Section 5.1. For EEWWG, 
which does not include any FSR, the expectation for the contribution from FSR from leptons 
(which is considered signal) is taken from PHOTOS. The experimental results correspond to a 
measurement with 10 % precision of the W + W~7 cross-section. The best agreement is obtained 
with KandY and RacoonWW; however, the measurements are of insufficient statistical precision 
to distinguish between the different calculations. The OPAL result is two standard deviations 
below the prediction of KORALW. Although the statistical significance is low, the 0(a) NL 
electroweak corrections of YFSWW implemented in KandY improve the agreement between 
data and Monte Carlo (the dominant effect is the inclusion of radiation from the W-bosons, 
specifically its interference with ISR). 





Data/Theory 


KandY 
RacoonWW 
EEWWG 
KORALW 


0.99 ±0.09 ±0.04 
0.98 ±0.09 ±0.06 
0.91 ±0.09 ±0.04 
0.84 ±0.08 ±0.04 



Table 3: The ratios of the experimental to expected SM W + W~7 cross-sections averaged 
over \fs for four theoretical calculations. The errors are from the statistical and systematic 
uncertainties on the measurement of the W + W~7 cross-section. 

New physics could appear as resonant structure in the W7 invariant mass distribution (for 
example the decay of an excited W-boson, W* — > W7). To investigate this possibility, for qq£z7f7 
and qqqq7 candidates the invariant masses of the two W ± 7 combinations in selected W + W~7 
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events are obtained from an additional kinematic fit. The fit uses the constraints of energy 
and momentum conservation and the constraint that the invariant masses of the reconstructed 
fif 2 and fsf 4 systems are both equal to the W-mass (previously the requirement was that 
both masses be equal). Only events for which the kinematic fit converges are retained. For MC 
events this cut rejects approximately 16 % of selected signal events. The W7 invariant mass 
is calculated from the four- momenta of the four fermions and the photon returned by the fit. 
Figure 4 shows the reconstructed invariant mass distribution for the two W ± 7 combinations 
for selected W + W~7 events with E 1 > 2.5 GeV. No resonant structure is observed. The data 
from the region |cos# 7 | < 0.80, where any contribution from new physics might be expected to 
be most apparent are also shown. 

5.1 Systematic Uncertainties 

The contributions to the systematic uncertainties on the W + W~7 cross-sections for the five 
values of y/s are listed in Table 4 and are described below. The total systematic errors are taken 
as the sum in quadrature of these components. When determining the average ratio of data to 
MC the systematic error components for the five energies are taken to be 100 % correlated. 



Systematic Uncertainty 


on (T WW7 /fb 






Error Source 


Variation 




< v/i >/GeV 








183 


189 


195 


201 


206 


Photons from jets 


±15 % 


9 


10 


12 


13 


13 


Photon energy scale 


±4 % 


6 


8 


5 


9 


10 


Photon angular acceptance 


±5 mrad 


4 


6 


4 


7 


7 


Photon energy resolution 


±10 % 


3 


4 


3 


5 


5 


W+W" Selection 


±1.1 % 


3 


4 


3 


5 


5 


Photon Identification 


±1.0 % 


3 


4 


3 


5 


5 


Photon Isolation 


±1.0 % 


3 


4 


3 


5 


5 


qq Background 


±6.5 % 


2 


2 


3 


3 


3 


Kinematic Fits 


±0.5 % 


1 


2 


1 


2 


2 


Monte Carlo Statistics 


±0.4 % 


1 


2 


1 


2 


2 


Luminosity 


±0.3 % 


1 


1 


1 


1 


1 


Total Systematic Error 




13 


17 


15 


20 


21 


Statistical Error 




117 


74 


84 


100 


73 



Table 4: The contributions to the experimental error on the W + W 7 cross-section for the five 
different values of y/s. The systematic variations on the various sources of error are indicated. 

Modelling of photons from jets: The modelling of photon candidates associated with the 
hadronic jets (both from FSR and from 7r° and rj decays) is studied by comparing the rate at 
which photons are identified in Z° — > qq events to the PYTHIA prediction (for this comparison 
data recorded at \fs ~ M z o during the 1998 — 2000 operation of the LEP accelerator are used). 
For 2.5 GeV< E 1 < 20 GeV, there are (38 ± 2) % more photon candidates identified in the 
data than expected from the Monte Carlo. Above 20 GeV the data are consistent with the 
Monte Carlo expectation. The ratio of data to Monte Carlo is used to estimate an energy- 
dependent correction (in photon energy bins of 2.5 GeV) to the Monte Carlo expectation for 
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the background from W + W~ events with photons associated with jets. After the W + W~7 
event selection, this corresponds to a (30 ± 2) % correction to the background from photons 
from jets 4 . Half the size of the correction is propagated as a systematic uncertainty. In the 
evaluation of the other systematic uncertainties all comparisons between data and MC are 
performed after making this correction. 

ECAL energy scale: A bias in the energy scale for photons (data relative to Monte Carlo) 
in the region of the energy cut, i.e. E 1 ~ 2.5 GeV, would result in a systematic bias in the 
W + W~7 cross-section measurement. The uncertainty on the ECAL energy scale for photons 
in this region is estimated by examining photons from n° decays in e + e~ — > qq events recorded 
at y/s ~ M z o during 1998—2000 and e + e~ — > qq(7) events recorded at y/s > 180 GeV. The 
mean reconstructed 7r° mass for 7r° candidates containing a photon with 2 GeV < E 1 < 3 GeV 
is (142 ± 2) MeV/c 2 in data compared to 137 MeV/c 2 in Monte Carlo. As a result a 4 % 
systematic uncertainty on the ECAL energy scale in the region of E 1 ~ 2.5 GeV is assigned. 
The resulting systematic uncertainty on the cross-section is 2 %. 

Photon Angular Acceptance: The systematic error associated with the requirement of 
|cos# 7 < 0.975| depends on the accuracy of the Monte Carlo simulation of the angular recon- 
struction from ECAL clusters at the edge of the acceptance. By comparing the reconstructed 
polar angle from different detectors (ECAL, tracking, muon chambers) the ECAL acceptance 
is known to ±3 mrad out to |cos# 7 | < 0.96. Beyond the tracking acceptance it is not possible 
to make this comparison. Therefore a 5 mrad uncertainty on the edge of the acceptance is 
assigned. 

As a cross-check a sample of ISR photons from e + e~ — >■ qq(7) events is used. Multi- 
hadronic events recorded at 180 GeV < y/s < 209 GeV are selected [25]. Photons are identified 
using the same criteria as for the W + W~7 cross-section analysis. In the data 241 photons are 
reconstructed in the region 0.950 < |cos 7 | < 0.975 compared to the Monte Carlo expectation of 
237.1. A 5 mrad bias between data and Monte Carlo would result in an expected discrepancy 
of 28.5 events in this region. The good agreement between data and Monte Carlo provides 
confirmation that the assigned uncertainty of 5 mrad is reasonable. 

ECAL energy resolution: The systematic error from the uncertainty in the ECAL energy 
resolution is obtained in a similar manner as that used for the ECAL energy scale using the 
same ir° sample. There is no evidence for a difference between data and Monte Carlo within 
the statistical precision of the comparison (±10 %). The precision of this comparison is used to 
assign a (10 %) uncertainty the energy resolution, which, when propagated to the uncertainty 
on the W + W~7 cross-section yields a systematic error of ±1 %. 

W+W - selection efficiency: Systematic uncertainties in the W + W~ event selection will 
result in corresponding uncertainties in the W + W~7 event selection. The estimated systematic 
uncertainty on the W + W~ selection efficiency is 1.1 % [1], where the largest uncertainties are 

4 For the comparison with Racoon WW given in Table 3 the systematic errors from photons from jets are 
calculated differently. In RacoonWW it is impossible to separate photons from FSR from quarks from other 
diagrams. Consequently the signal definition is modified to include all FSR photons within the theoretical 
acceptance cuts. In this case the data/MC discrepancy for photons from jets in Z° — > qq may either be assigned 
to a mis-modelling of FSR (signal) or to a mis-modelling of hadron production rate (background) . Consequently 
the systematic uncertainties are larger than for the case when FSR from quarks is also treated as background. 
The central value for the RacoonWW comparison uses the average of the results obtained and half the difference 
is assigned as a systematic error. 
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related to the QCD and fragmentation modelling of jets. For the data sample considered here, 
the W + W~ event selection yields 11752 events which is statistically compatible with the Monte 
Carlo expectation of 11670 ± 58 (where the error is taken to be the theoretical uncertainty on 
the CC03 cross-section). The difference is consistent with the quoted systematic error of 1.1 %. 

Photon Identification : A systematic uncertainty of 1 % is assigned to cover the uncertain- 
ties in the simulation of the photon conversion rate and the accuracy of the simulation of the 
electromagnetic cluster shape [26]. Systematic uncertainties arising from the isolation require- 
ments are discussed below. The efficiency obtained from KandY is consistent with that from 
KORALW and no additional systematic uncertainty is assigned. 

Photon Isolation : The systematic error associated with the isolation requirements depends 
on the accuracy of the Monte Carlo simulation of the fragmentation process in hadronic jets. 
This is verified in Z° — > qq events recorded at y/s ~ M z o during 1998 — 2000. For each selected 
event, the inefficiency of the isolation requirements is determined for random orientations of 
the isolation cone and parametrised as a function of the angle between the cone and the nearest 
jet. For all jet-cone angles the inefficiency in the Monte Carlo and data agree to better than 
1 %, consequently a 1 % systematic error is assigned. Consistent results, albeit with lower 
statistical precision, are obtained from W + W~ — ► qqiV e events. 

As a cross-check of the photon identification and isolation requirements the sample of recon- 
structed photons in e + e~ — > qq(7) events is used. The ratio of the number of reconstructed pho- 
tons with 2.5 GeV < E 1 < 50 GeV in the data to the Monte Carlo expectation is 1.015 ±0.023. 
Good agreement is observed over all cos# 7 . Due to the limited statistical sensitivity of this test 
no additional systematic uncertainty is assigned to the photon identification/isolation efficiency 

qq~7 background: The dominant source of non-W + W~ background is from e + e~ — > Z°/7 — > 
qqy where the identified photon candidate is a genuine photon from 1SR. Uncertainties in the 
modelling of QCD/fragmentation lead to systematic uncertainties in the level of background 
from e + e~ — ► qqy events in the W + W~ event selection [1]. As a result the qqy background 
in the W + W~y selection is uncertain to 6.0 %. An additional systematic error of 2.5 % arises 
from the uncertainties in the modelling of ISR in e + e~ — > qqy events. 

Kinematic Fits : The W + W~y event selections require that a kinematic fit converges and 
has a reasonable probability. Possible mis-modelling of the detector response/resolution could 
result in a difference in the rates at which the fits fail for data and Monte Carlo. This was 
checked by applying the kinematic fits used in the W mass analysis to all selected W + W~ 
events and comparing the failure rates for data and Monte Carlo. The efficiency (qqiVi and 
qqqq combined) in data is (83.9 ±0.4) % compared with the Monte Carlo expectation of 84.1 %. 
The ratio of these efficiencies is 0.997 ± 0.005 and, consequently, a systematic uncertainty of 
0.5 % is assigned. 

Monte Carlo Statistics : The effect of finite Monte Carlo statistics is taken into account and 
leads to 0.3 % systematic uncertainties on the measured cross-sections. 

Luminosity: The total uncertainty on the integrated luminosity of the data samples is 0.3 %, 
dominated by systematics. 
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6 Limits on Mw ®{pt) Systematic Uncertainties 



The anticipated experimental error on Mw from LEP2 is approximately 35 MeV. A potential 
source of theoretical uncertainty is the treatment of higher order QED corrections in the Monte 
Carlo programs used to simulate the process e + e~ — > 4/(7). A recent estimate suggests a 
total theoretical systematic uncertainty due to 0(a) effects of 5 MeV [7]. However, as pointed 
out by the authors [7] , this estimate is based upon the invariant mass of the yTV^ system in 
e + e~ — > /i~V IJi ud('j) events, whereas the experimental procedure used to extract Mw is compli- 
cated by the fact that the four LEP collaborations use kinematic fits to improve significantly 
the event-by-event W-mass resolution [27-30]. One effect of the kinematic fit is to constrain 
the total energy of the reconstructed fermions to ^/s. For events with photons from ISR this 
procedure introduces a bias in the reconstructed W-mass as the energies of four fermions should 
be constrained to V^, the centre-of-mass energy after photon radiation, rather than to y/s. 
Consequently, as a result of the experimental procedure used to extract Mw, the 0(a) theoret- 
ical systematic uncertainties may be significantly greater than those obtained by considering 
the invariant mass distribution of the final state fermions [31]. 

6.1 QED and Electroweak Corrections in KandY 

In the KandY generator it is possible to study the effects of different theoretical corrections 
using event correction weights [9] which, when used to weight generated events, allow different 
theoretical predictions to be tested. By processing generated fully-simulated events through 
the full OPAL W-mass analysis it is possible to determine the W-mass biases associated with 
these corrections. For example, degrading the 0(a 3 ) exponentiated LL treatment of collinear 
ISR to 0(a 2 ) results in a systematic bias of less than 1 MeV [7,30]. In a similar manner the 
non-leading (NL) 0(a) electroweak corrections, including radiation from the W-bosons, may 
be switched off using the appropriate event correction weights, Wj^- When applied to the full 
OPAL W-mass analysis it is found that dropping the O(a) NL electroweak corrections results 
in a shift in the reconstructed W-mass of 15 MeV. This relatively large bias is due to the 
modification of the \fs i distribution rather than a distortion in the invariant mass distribution 
of the fermion pairs [31]. The change in the yfs' distribution is due to the inclusion in KandY 
of the diagrams for radiation from the W-bosons which, through interference with the ISR 
diagrams, reduces the cross-section for the production of real photons [32] 5 . Although the 
fractional change in W + W~7 cross-section is largest at cos# 7 = where the photon production 
rate is reduced by 30 % [9], in absolute terms the reduction in the cross-section shows no strong 
cos# 7 dependence. Consequently the W + W~7 cross-section measurement provides a test of the 
modelling of radiation from the W-bosons (and the interference with ISR) in the KandY Monte 
Carlo. Since the largest source of systematic bias from the so-called 0(a) NL corrections is 
a direct result of the modification of the spectrum of real photons, the associated systematic 
uncertainties on M w may be constrained by the measurement of the W + W~7 cross-section. 

6.2 Constraints from the W + W~7 Measurements 

To investigate the experimental limits on possible biases on the measurement of M w due to 
photon production away from the collinear region the correction weights from KandY are first 

5 This effect was investigated by running YFSWW with KEYC0R=2 and KEYCDR=3 switching between the YFS 
form factor solely for ISR and the full form factor including WSR and interference between WSR and ISR. 
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expressed in the form 

By modifying the weights to 

it is possible to investigate a continuous range of scenarios. A value of k = corresponds to 
the treatment of 0(a) NL electroweak corrections of YFSWW {i.e. the default in KandY) and 
k = 1.0 corresponds to dropping the NL 0(a) corrections. The parameter k and its errors are 
obtained from a binned extended maximum likelihood fit to the |cos# 7 | distribution of Figure 2, 
taking into account both the overall normalisation and shape, giving 

k = 0.38 ±0.45 ±0.15, 

where the first error is statistical and the second due to systematic uncertainties in the event 
selection efficiency The data favour the KandY prediction including the NL corrections. Most 
of the sensitivity comes from the photon rate rather than the angular distribution. The mea- 
sured value of k suggests that the measured value of M w from the OPAL W-mass analysis, 
obtained using KandY as a reference, should be corrected by (—5 ±6) MeV. Using the measured 
cross-section alone gives a similar result of (— 1 ± 7) MeV. From these studies it is concluded 
that the systematic error on Mw due to the Monte Carlo implementation of QED diagrams 
resulting in real photon production away from the collinear region should be not more than 

6 MeV. 

7 Anomalous Quartic Gauge Boson Couplings 

The non-Abelian nature of the electroweak sector of the Standard Model results in vector 
boson self-interactions. In addition to the triple gauge boson couplings (TGCs), W + W~7 
and W + W~Z°, the Standard Model predicts the existence of four quartic gauge couplings, 
W+W-W+W", W+W-Z°Z°, W+W-Z°7 and W+W"77. These couplings are not expected to 
play a significant role at LEP energies, but will be important at the LHC [33] and at a future 
TeV linear collider [34]. 

Quartic gauge boson couplings can be probed in final states with three vector bosons. At 
LEP centre-of-mass energies, final states involving three massive gauge bosons are kinematically 
out of reach. However, it is possible to study the processes e + e~ — > W + W~7 [10, 11] and 
e + e~ — > Z°77 [35]. In the Standard Model, the contribution of the quartic couplings to e + e _ — > 
W + W~7, shown in Figure 5, is expected to be too small to measure and that to e + e" — > Z°77 
is zero. Nevertheless, it is possible to set direct limits on possible anomalous contributions to 
the quartic gauge boson couplings. 

7.1 Theoretical Framework 

In the SM the form and strength of the vector boson self-interactions are fixed by SU(2)xU(l) 
gauge invariance. As is the case for triple gauge boson couplings [36], in extensions to the SM, 
anomalous quartic couplings can be parametrised by additional terms in the Lagrangian [20, 
37,38]. These are required to conserve custodial SU(2) C symmetry in order to avoid deviations 
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of the p parameter 6 from the experimentally well established value close to 1. Only operators 
which do not introduce anomalous triple gauge couplings are considered. For example, the 
anomalous quadrupole moment operator generates both W + W~7 and W + W~77 couplings. 
Therefore, it is not considered as a source of genuine anomalous quartic couplings since its 
strength, A 7 , is already tightly constrained from the study of TGCs at LEP [39,40] and at 
the Tevatron [41]. The lowest dimension operators which generate genuine anomalous quartic 
couplings involving photons are of dimension six. Three such possibilities are considered here, 
£°, C% [37] and Q [20,42]: 

^6 

^6 
rn 

with 

V ZJ cos 9 W J 

where F Miy and W^ u are the field strength tensors of the photon and W fields respectively. 
Both £° and >Q, which conserve C and P (separately), generate anomalous W + W~77 and 
Z°Z°77 couplings. The CP violating term £g results in an anomalous W + W~Z°7 coupling. 
In each case, the strength of the coupling is proportional to a^/A 2 , where A represents a scale 
for new physics. A more general description of the operators leading to anomalous quartic 
couplings accessible at LEP can be found in the paper of Belanger et al. [43] . The two additional 
dimension 6 operators, parametrised by ao and a c , identified by Denner et al. [44] are not 
considered here as the effects of ao and a c are almost identical to those of ao and a c , respectively. 

7.2 Experimental Limits 

The selected W + W~7 events are used to set limits on possible anomalous contributions to the 
W+W"77 and W + W"Z°7 quartic gauge couplings. The limits are extracted from the measured 
differential cross-section as a function of the photon energy and photon polar angle. The signal 
of anomalous quartic gauge boson couplings at LEP would be an excess of W + W~7 events. 
The effect of anomalous QGCs increases with photon energy. Furthermore, the sensitivity to 
anomalous QGCs increases with increasing y/s. 

The calculation of Stirling and Werthenbach [20] allows for the assessment of the impact of 
anomalous quartic couplings and is implemented in the EEWWG program. This calculation 
includes the ISR diagrams, the WSR diagrams, the SM QGC diagram and can accommodate 
anomalous quartic couplings. However, the recent implementation of anomalous QGCs in the 
Racoon WW [44] and WRAP [45] programs identified a problem with the EEWWG program, 
indicating that a — > — a and a c — > — a c in EEWWG. In this study the EEWWG program is 
used with the signs of a and a c inverted. To set limits on possible anomalous couplings a binned 

6 p = M^ V /(M^„ cos 2 9 W ), where M w and M z o are the masses of the W^ 1 and Z° bosons and 9w is the weak 
mixing angle. 
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maximum likelihood fit to the observed distribution of [E 7 , |cos# 7 |] is performed using bins of 
[5 GeV, 0.1]. Fits are performed to the data for the five separate energy ranges of Table 1 and 
the resulting likelihood curves are summed. The effects of anomalous couplings are introduced 
by reweighting events generated with KandY using the average ratio of anomalous QGC to 
SM matrix elements from EEWWG in the relevant bin of [E 1 , |cos# 7 |]. The resulting summed 
likelihood curves are shown in Figure 6. Results are obtained for three single parameter fits, 
where one of do, a c or a n is varied whilst the other two parameters are set to zero, and a two 
parameter fit to {a ,a c }. The results include the effect of the experimental systematic errors 
and assume a 10 % theoretical uncertainty 7 on the cross-section for e + e~ — > W + W"7. These 
uncertainties are taken to be 100 % correlated between the five energy ranges. The best fit 
does not occur at the SM value of zero. However this does not imply the data are inconsistent 
with the SM. The consistency with the SM prediction, given by the probability of obtaining 
a value of — ln£ greater than that observed for {ao = 0,a c = 0,a n = 0}, is 19 %. The 95 % 
confidence level upper limits on the anomalous couplings, obtained from the likelihood curves, 
A(ln£) = 1.92, are: 

-0.020 GeV" 2 < a /A 2 < 0.020 GeV" 2 , 
-0.053 GeV" 2 < a c /A 2 < 0.037 GeV" 2 , 
-0.16 GeV" 2 < a n /A 2 < 0.15 GeV" 2 , 

For a c the region -0.020 GeV" 2 < a c /A 2 < -0.002 GeV" 2 is also excluded at the 95 % C.L. 
The derived upper limits are less restrictive than the expected limits. For example, the expected 
limit on a is |a /A 2 | < 0.014 GeV" 2 . The limits are worse than expected due to a slight excess 
of high energy photons in the y/s > 205 GeV data sample. 

8 Conclusions 

Using 187 W + W~7 candidates with photon energies greater than 2.5 GeV the W + W~7 cross- 
section is measured at five values of yfs. The results are consistent with the Standard Model 
expectation. Averaging over the five energies, the ratio of the observed cross-section to the 
prediction of the concurrent Monte Carlo KoralW and YFSWW (KandY) is 

^(data/MC) = 0.99 ± 0.09 ± 0.04, 

where the errors represent the statistical and systematic uncertainties respectively. This pro- 
vides a 10 % test of the KandY implementation of 0(a) effects producing a real photon away 
from the collinear region. From these studies it is concluded that the systematic error on M w 
due to the Monte Carlo implementation of QED diagrams resulting in real photon production 
away from the collinear region should be not more than 6 MeV. 

The data are used to derive 95 % confidence level upper limits on possible anomalous 
contributions to the W + W~77 and W + W"Z°7 vertices: 

-0.020 GeV" 2 < a /A 2 < 0.020 GeV" 2 , 
-0.053 GeV" 2 < a c /A 2 < 0.037 GeV" 2 , 
-0.16 GeV" 2 < a n /A 2 < 0.15 GeV" 2 , 

where A represents the energy scale for new physics. 

7 This represents a conservative estimate of the theoretical uncertainty, comparisons of YFSWW and 
RacoonWW suggest 5 %. 
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10 Appendix 

For the purpose of the combination of results from the four LEP experiments cross-section 
results are obtained for the signal definition: 

• E~ f > 5 GeV, 

• |cos# 7 | < 0.95, 

• cos# 7/ < 0.90, 

• \M f j 2 - M w | and \M f .j 4 - M w | < 2T W . 

The experimental cuts on the photon acceptance are modified to match the signal definition. 
For the modified selection, 124 events are selected, compared to the SM expectation (KandY) 
of 118.7 ± 0.6. The results are summarised in Table 5. 
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o WW7 / fb 


< v/i >/GeV 


Data 


KandY 


182.68 


102 ± 60 ±5 


141 ±2 


188.63 


163 ± 41 ±6 


175 ±3 


194.44 


166 ± 57 ±7 


201 ±2 


200.21 


214 ± 60 ±7 


216 ±3 


205.92 


298 ± 50 ±8 


226 ±3 



Table 5: W + W~7 cross-section measurements for the signal definition to be used for a LEP 
combination of results. The errors are statistical and systematic respectively. The systematic 
uncertainties are calculated as described in Section 5.1. 
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OPAL 




Figure 1: The five kinematic variables used to classify the photon in W + W~7 events as being 
from ISR, FSR or associated with the jet. Unless otherwise specified the distributions are shown 
for qq£z7^ and qqqq events combined. The variables are: a) the angle between the photon and 
the nearest jet, cos# 7 _ JET ; b) the angle between the photon and the charged lepton, cos 6^ 
(qq£l^ events only); c) the reconstructed W-boson mass under the hypothesis that the photon is 
associated with jet; d) the reconstructed W-boson mass under the hypothesis that the photon is 
from ISR; e) the reconstructed W-boson mass under the hypothesis that the photon is from FSR 
(qq£z7^ only). Plot f) shows the resulting relative likelihood distribution for the ISR hypothesis. 
In all cases the data are shown by points with error bars, the total SM expectation is shown by 
the histogram and the contributions from processes other than ISR are shown by the hatched 
histograms. 
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Figure 2: For selected W + W"7 events (180 GeV < y/s < 209 GeV), a) shows the photon 
energy spectrum, b) the modulus of the cosine of the polar angle of the photon, and c) the 
cosine of the angle between the photon and the nearest charged fermion. The data are shown 
by the points with error bars and the SM expectations (KandY) are shown by the histograms. 
The doubly-hatched histograms indicate the contributions from non-W + W~ background and 
background from photons associated with the parton-shower (either FSR or from hadron decay). 
The singly-hatched histograms show the contributions from FSR from leptons. The expected 
E 1 and |cos# 7 | distributions for an anomalous QGC of a /A 2 = 0.040 GeV -2 are also shown. 
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Figure 3: Measured W + W~7 cross-section for the signal definition of Section 3. The points 
with error bars show the OPAL measurements. The curve shows the SM expectation obtained 
from the KandY program. 
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Figure 4: Reconstructed invariant mass of W ± 7 in selected W + W~7 events with > 2.5 GeV 
(two entries per event). The data are shown by the points, the Standard Model expectation, 
determined from KandY, is shown by the histogram. The singly hatched histograms show the 
contribution from FSR from leptons and the doubly hatched histograms show the background. 




Figure 5: Standard Model production diagram for the W + W 7 final states involving the 
W + W~77 and W + W~Z°7 quartic gauge couplings. 
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Figure 6: Likelihood curves for the anomalous QGC parameters a , a c and a n . Also shown is 
the 95 % C.L. region for (ao,a c ). The curves include the experimental systematic uncertainties 
and a 10 % theoretical uncertainty for the e + e~ — > W + W~7 cross-section. 
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